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ABSTRACT 
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Reflected  shock  waves  were  used  to  heat  a variety  of  formaldehyde  mixtures  to — 

18  3 

temperatures  of  1600-3000K  at  total  concentrations  near  5x10  /cm  . Formaldehyde 
decay  was  monitored  by  IR  emission  at  3.50  ym  using  both  O2  and  ^0  as  the  oxidant 
In  mixtures  highly  diluted  with  argon.  Other  mixtures  with  added  CO  were  monitored 
for  oxygen  atom  production  via  the  flame-band  emission  at  450  nm  and  for  CO2  produc- 
tion via  IR  emission  at  4.27  ym.  These  data  were  then  compared  to  the  results  ob- 
tained from  numerical  integration  of  the  rate  equations  based  upon  a likely  mechan- 
ism. It  was  possible  to  achieve  good  agreement  with  the  formaldehyde  reaction  rate 
constants  listed  below: 


(1)  CH20  + M = HCO  + H + H,  k1  = 5.5xl0"8  exp(-339  kJ/RT)  cm3  s-1 

(2)  CH20  + H = HCO  + H2,  k2  = 5.5xl0"10  exp(-43.9  kJ/RT)  cm3  s_1 

(3)  CH20  + OH  = HCO  + H20,  k3  = 1 .25xl0-1 1 exp(-0.7  kJ/RT)  cm3  s"1 

(4)  CH20  +0  = HCO  + OH,  k4  = 3.0x10  exp(-12.9  kJ/RT)  cm3  s"1 

The  values  listed  for  k^  and  k4  are  taken  directly  from  the  low  temperature  studies 
of  these  reactions;  the  k^  and  k2  values  also  fit  recent  experiments  in  this  labora- 
tory on  CH2O  pyrolysis.  The  ^0  experiments  also  served  as  a check  on  the  rate  con- 
stant assignment  for  the  reaction: 

(19)  H + N20  = N2  + OH,  klg  = 1.5xl0-9  exp(-92  kJ/RT)  cm3  s"1. 

Reasonable  agreement  was  achieved  here.  Sensitivity  studies  would  suggest  these 
rate  constants  are  probably  accurate  to  within  a factor  of  two  between  1700-2500K. 

The  values  of  kj  and  k^  obtained  in  this  work  are  at  least  an  order  of  magni- 
tude lower  than  the  values  frequently  used  in  current  modeling  studies  of  methane 
oxidation.  The  implications  of  these  differences  are  discussed  in  light  of  recent 
suggestions  that  C2  species  participate  In  methane  oxidation. 
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INTRODUCTION 


Formaldehyde  has  long  been  recognized  as  an  important  intermediate  in  many 
hydrocarbon  oxidation  systems.  As  a result,  it  is  usually  included  in  the  various 
oxidation  models.  Unfortunately,  remarkably  little  direct  high  temperature  data  is 
available  on  this  system;  the  rate  constant  assignments  are  therefore  uncertain. 
This  in  turn  introduces  a considerable  ambiguity  into  mechanistic  analysis  of  even 
simple  hydrocarbon  oxidation  systems. 

To  amplify  these  comments,  consider  that  the  recent  review  of  Jensen  and  Jones 
[1]  suggests  that  the  reactions 


CH20  + OH  - HCO  + H20 
CH20  + H -*■  HCO  + H2 
CH20  +0  - HCO  + OH 

have  uncertainty  factors  of  100,  30,  and  30,  respectively.  (Here  the  uncertainty 
factor  UF  is  defined  such  that  k x UF  and  k/UF  provide  probable  approximate  upper 
and  lower  bounds  to  the  rate  constant  k.)  Similarly,  shock  tube  and  flame  studies 
report  rate  constants  for  formaldehyde  dissociation  that  differ  by  several  orders 
of  magnitude  [2,3].  In  this  light,  it  is  not  surprising  that  some  combustion  mod- 
elers have  advocated  that  additional  experiments  need  to  be  done  on  the  CH20  sys- 
tem [4]. 

To  date,  there  have  been  only  three  shock  tube  studies  which  have  focused  di- 
rectly upon  the  formaldehyde  reactions.  The  pioneering  study  of  both  pyrolysis  and 
oxidation  was  done  by  Gay,  et  al.  [5].  Several  chain  mechanisms  with  hydrogen  atom 
Intermediates  were  proposed  to  explain  the  pyrolysis  results,  but  all  predicted  ac- 
tivatlon  energies  significantly  higher  than  observed.  Thus  it  was  not  possible  to 
assign  any  elementary  reaction  rate  constants.  The  oxidation  studies  suggested  that 
formaldehyde  decomposition  is  greatly  accelerated  by  the  addition  of  oxygen.  Schecker 
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and  Jost  [2]  studied  the  pyrolysis  reaction  between  1400-2200K  and  varied  the 
CH2O  mole  fraction  from  0.02 % to  1%.  The  data  were  reasonably  consistent  with 
the  earlier  study.  They  were  analyzed  in  terms  of  a chain  mechanism,  and  a rate 
constant  was  reported  for  the  dissociation  step.  The  one  troubling  aspect  of  this 
analysis  was  that  the  activation  energy  for  the  1%  mixtures  was  lower  than  one  would 
expect  from  the  steady  state  analysis.  Recently,  studies  in  this  laboratory  of  the 
pyrolysis  reactions  used  an  extensive  numerical  integration  study  as  part  of  the 
analysis  [6],  Here  it  was  not  necessary  to  invoke  the  steady  state  assumption. 

It  was  shown  that  the  mechanism  used  by  Schecker  and  Jost  was  consistent  with  the 
observations,  but  that  the  rate  constant  assignments  must  be  modified.  The  decay 
data  were  influenced  primarily  by  the  CH^O  dissociation  reaction  and  by  the  hydro- 
gen atom  attack  upon  (^O.  It  was  not  possible  to  completely  separate  the  effects 
of  these  reactions,  and  thus  unambiguous  assignments  were  not  possible.  However,  it 
was  evident  that  the  dissociation  rate  constant  was  much  lower  than  that  reported 
by  Schecker  and  Jost. 

The  present  CH^O  oxidation  study  was  initiated  for  two  reasons:  (1)  To  attempt 
to  complement  the  pyrolysis  study  by  separating  the  contributions  from  the  reactions 
CH^O  + M -*•  HCO  + H + M and  CHgO  + H -*■  HCO  + H2»  thus  allowing  rate  constant  assign- 
ments to  be  made.  (2)  To  obtain  more  direct  high  temperature  information  about  the 
reactions  CHgO  + 0 -►  HCO  + OH  and  ^0  + OH  -*■  HCO  + H2O.  The  approach  utilized  is 
similar  to  that  used  previously;  concentration-time  behavior  is  monitored  and  com- 
pared to  that  calculated  by  numerical  integration  of  the  assumed  mechanism.  Several 
techniques  were  employed  to  attempt  to  minimize  ambiguity  in  rate  constant  assign- 
ments: (1)  Concentration-time  profiles  over  large  extents  of  reaction  were  ob- 
tained for  three  species  --  CHgO,  CO2,  and  0.  (2)  Conditions  were  chosen  to  mini- 

mize the  number  of  unknown  variables.  (3)  Additional  experiments  with  ^0  as  oxi- 
dant were  used  to  confirm  assignments  made  with  02  as  oxidant;  the  radical  concen- 
trations differ  significantly  in  the  two  systems. 


With  this  approach  It  was  possible  to  analyze  the  data  to  yield  rate  constant 
assignments  for  all  four  CH20  reactions.  The  N20  system  also  allowed  a check  on 
the  rate  constant  for  H + N20  -*■  OH  + N2<  Perhaps  the  most  significant  conclusion 
of  this  work  is  that  several  of  the  high  temperature  formaldehyde  reaction  rate  con- 
stants appear  to  be  much  lower  than  those  currently  utilized  in  modeling  calculations 


EXPERIMENTAL 


The  7.6  cm  shock  tube,  gas  handling  system,  and  optical  configuration  have 
been  described  previously  [11].  Infrared  emissions  were  collimated  by  two  slits 
1.5  mm  wide  and  5 mm  high  placed  50  mm  apart.  After  passing  through  an  appropriate 
filter,  the  emissions  were  focused  upon  a liquid-nitrogen-cooled  indium  antimonide 
detector;  the  detector-preamp  combination  had  a nominal  rise  time  of  1 ps.  To  fol- 
low CH2Q,  either  a 3.50  pm  interference  filter  (FWHM  = 0.12  ym)  or  a filter  wheel 
set  at  3.50  ym  (FWHM  = 0.15  ym)  was  used.  C02  was  monitored  with  a 4.27  ym  inter- 
ference filter  (FWHM  = 0.18  ym).  The  flame-band  emissions  were  monitored  by  a RCA 
1P28A/V1  photomultiplier  mounted  behind  a 450  nm  interference  filter  (FWHM  = 6.5  nm) 
and  two  slits  1 mm  wide  by  2 mm  high  which  were  38  mm  apart.  The  measured  response 
time  of  the  photomultiplier/preamp  system  was  less  than  2 ys.  Data  were  collected 
with  either  a Biomation  Model  805  transient  recorder  or  the  system  previously  de- 
scribed [12].  Biomation  data  was  typically  collected  at  0.5  ys  intervals  while  the 
other  device  collected  data  at  10.0  ys  intervals. 

Formaldehye  was  prepared  by  heating  paraformaldehyde  to  ^ 100°C,  passing  the 
distillate  through  a dry  ice-acetone  trap,  and  then  condensing  the  monomer  at  77K. 
Care  was  taken  to  keep  the  monomer  pressure  below  1 kPa  (7.5  torr)  at  all  times  to 
minimize  polymerization.  The  monomer  was  then  subjected  to  several  bulb-to-bulb 
distillations  from  dry  ice  to  liquid  nitrogen  temperatures.  Then,  while  at  -78°C, 

It  was  distilled  Into  the  mixing  bulb.  A fresh  batch  of  CH20  was  prepared  for  each 
mixture.  Other  gases  used  were  AIRCO  Ar  (99.9995%)  and  N20  (99.995%),  Matheson  CO 
(99.99%),  and  Scientific  Gas  Products  02  (99.999%).  The  CO  was  slowly  passed  through 
a coll  of  copper  tubing  at  77K  before  addition  to  the  vacuum  line.  Other  gases 
were  used  as  supplied. 


The  shock  tube  was  pumped  down  to  'v  3 mPa  and  the  observed  leak-outgassing  rate 


was  usually  near  4 mPa/min,  The  tube  was  isolated  from  the  pumps  for  approxi- 
mately one  minute  prior  to  shock  initiation;  the  nominal  background  pressure  was 
near  7 mPa.  Since  the  test  section  was  pressurized  to  2.7  kPa  with  the  mixtures 
of  interest,  the  background  impurities  from  the  shock  tube  were  present  in  the  low 
ppm  range  in  the  shocked  mixtures.  Mylar  diaphragms  and  helium  driver  gas  were 
used  throughout. 

Temperature  behind  the  reflected  shocks  were  computed  in  the  usual  way  from 
measured  incident  shock  velocities.  Reflected  shock  pressures  were  measured  with 
a fast  response  pressure  transducer  and  were  always  found  to  be  in  good  agreement 
with  those  calculated.  Furthermore,  the  reflected  shock  pressure  was  always  con- 
stant over  the  time  interval  that  data  were  collected;  as  a result,  no  attempt  was 
made  to  correct  for  non-ideal  effects. 
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RESULTS 


Formaldehyde  Decay 

The  mixtures  studied  are  listed  in  Table  1.  All  data  were  collected  at  0.2- 
0.5  ys  intervals.  The  interference  filter  was  used  for  Mixture  A while  the  filter 
wheel  was  used  in  B.  The  light  transmitted  by  the  filter  was  approximately  six  times 
that  of  the  wheel;  the  signals  observed  in  A were  significantly  greater  in  spite  of 
the  lower  C^O  concentration.  In  Mixture  B the  emission  was  observed  to  decay  com- 
pletely, whereas  in  A the  signal  was  seen  to  decay  to  a small, approximately  constant 
signal  which  was  5 - 10%  of  the  peak  signal  observed.  The  background  emission  is 
probably  due  to  i^O  and/or  OH  [13].  No  such  signal  would  be  expected  for  B;  here 
the  lower  sensitivity  would  have  resulted  in  a background  signal  comparable  to  the 
noise  level.  Attempts  were  made  to  correct  the  signal  in  Mixture  A by  subtracting 
the  estimated  background  signal  from  the  total  observed.  For  simplicity,  it  was  as- 
sumed that  the  background  signal  rose  linearly  from  zero  at  t = 0 to  the  constant 
level  seen  shortly  after  the  total  signal  ceased  to  decay.  Results  utilizing  such 
a correction  are  shown  in  Fig.  1.  Fortunately,  the  shape  of  the  (^0  signal  is  such 
that  much  of  it  is  quite  insensitive  to  the  shape  assumed  for  the  background  emission 
For  example,  assuming  the  background  signal  to  be  a step  function  at  t = 0 yields 
a signal  virtually  indistinguishable  from  that  shown  in  Fig.  1 for  10  < t 5 50  ys; 
even  at  60  ys  this  other  version  of  the  correction  is  only  10%  lower.  In  spite  of 
this  relative  insensitivity,  it  was  felt  that  attempts  to  characterize  the  signal 
in  terms  of  an  exponential  decay  parameter  were  unwarranted;  not  only  could  small 
errors  in  the  correction  markedly  influence  the  apparent  decay  rate  since  this  is 
the  slope  measurement,  but  also  the  signal  only  appeared  to  approach  exponential 
decay  in  the  region  between  1.0  and  0.1  volts  where  there  are  obvious  S/N  problems. 

As  a result,  three  points  were  chosen  in  an  attempt  to  characterize  the  decay  profile 
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The  first,  the  induction  time,  t . , is  only  applicable  for  the  lower  temperature  ex- 
periments where  there  was  an  obvious  plateau  in  the  signal  prior  to  the  onset  of 
decay.  Here  t^  was  defined  as  the  time  at  which  the  corrected  signal  had  decayed 
5%  from  its  plateau  value.  For  convenience  in  the  data  reduction,  these  times  are 

measured  relative  to  t , the  time  at  which  the  reflected  shock  was  first  visible  to 

o 

the  detector.  This  point  is  typically  5 ys  earlier  than  the  time  of  shock  passage 
of  the  window  mid-point.  (The  consequences  of  this  finite  slit  width  are  discussed 
later.)  At  higher  temperatures,  the  plateau  became  a peak  and  no  induction  times 
were  recorded.  Operationally,  the  signal  was  considered  to  exhibit  a plateau  if  it 
maintained  a constant  level  for  at  least  five  microseconds.  The  second  parameter 
used  was  t',  the  time  (again  relative  to  tQ)  at  which  the  corrected  signal  decayed 
to  3.0  volts.  Typically  this  parameter  characterizes  a point  early  in  the  decay 
profile  where  any  uncertainties  in  the  level  of  the  background  emission  should  have 
a negligible  effect.  The  final  parameter  is  t",  the  time  relative  to  t when  the 
corrected  signal  decayed  to  0.3  volts.  Use  of  this  parameter  allowed  for  measurement 
of  the  time  interval  corresponding  to  a decade  of  decay,  but  at  the  price  of  some  ad- 
ditional uncertainty  since  the  total  signal  at  this  point  was  only  two  to  three  times 
the  background  level.  These  three  parameters  are  listed  in  Table  1.  An  error  analy- 
sis of  these  parameters  considered  three  likely  sources  of  error:  (1)  Measurement 
errors  due  to  noise  in  a single  experiment  (including  additional  uncertainties  in 
t"  introduced  by  the  background  corrections).  (2)  Uncertainty  in  location  of  tQ: 

This  value  was  calculated  from  measured  incident  shock  velocities  to  avoid  bias; 
extensive  comparisons  of  calculated  values  versus  the  first  measurable  baseline 
deviations  when  step  function  signals  were  expected  (i.e.,  CO2  calibration  experi- 
ments) suggest  t was  known  to  ±2  ys.  (3)  Temperature  uncertainty:  Using  the  pro- 
cedures described  earlier  [11],  each  experiment  was  characterized  in  terms  of  both 
a constant  incident  shock  velocity  temperature  as  well  as  one  obtained  assuming  the 
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small  variations  in  time  intervals  between  various  velocity  stations  indicated  an 
accelerating  shock  front.  The  difference  in  these  two  temperatures  was  used  as  a 
measure  of  the  temperature  uncertainty  and  typically  ranged  from  zero  to  50K.  These 
differences  for  a given  mixture  were  averaged  to  give  an  overall  uncertainty.  The 
average  values  were  then  used  in  conjunction  with  "local  slope"  measurements  of  the 
various  parameter-temperature  plots  to  estimate  the  errors  likely  to  be  caused  in 
the  parameter  by  the  small  variations  in  temperature  expected.  The  three  likely 
sources  of  error  were  then  combined  in  a standard  propagation-of-error  treatment 
to  obtain  an  estimate  of  the  overall  error.  Plots  of  the  various  parameters  as  a 
function  of  temperature  yielded  curves  where  there  was  good  agreement  between  the  ob- 
served deviations  from  a fitted  line  and  the  error  bars  calculated  as  described 
above.  Thus  it  is  felt  the  parameters  chosen  are  reliable  indices  of  C^O  decay 
in  Mixture  A. 

Similar  parameters  were  used  to  characterize  Mixture  B.  The  only  differences 
here  were  choices  of  1.0  volt  and  0.1  volt  to  define  t1  and  t".  As  mentioned  earlier, 
the  smaller  observed  intensity  here  required  this  change.  There  was  no  background 
emission  problem  here,  and  semilog  plots  of  the  emission  signal  versus  time  often 
were  linear,  suggesting  exponential  decay.  However,  the  lower  region  of  the  plots 
invariably  had  an  appreciable  noise  component,  and  there  could  have  easily  been  cur- 
vature within  the  scatter.  As  a result,  it  was  felt  that  the  time  parameter  approach 
was  a better  method  to  attempt  to  characterize  the  data.  The  parameters  are  listed 
in  Table  1.  The  error  analysis  here  was  similar  to  that  employed  for  Mixture  A. 

Again  plots  of  the  parameters  versus  temperature  indicated  a degree  of  scatter  con- 
sistent with  the  calculated  error  bars.  It  should  be  noted  that  all  the  times  here, 

as  in  Mixture  A,  are  relative  to  t . 

o 

Note  in  Table  1 that  the  induction  times  are  appreciably  larger  in  the  Og  system 
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even  though  the  decay  rate,  once  reaction  is  initiated,  appears  to  be  comparable 
in  the  two  systems. 

Flame-band  Signal 

Information  about  oxygen  atom  production  was  obtained  by  observing  emissions 
at  450  nm  in  systems  to  which  CO  has  been  added.  The  procedure  is  completely  analo- 
gous to  that  used  previously  [14].  Mixtures  studied  are  listed  in  Table  2.  As  in 
earlier  work,  additional  experiments  were  performed  which  were  similar  to  those  in 
Table  2 except  that  CO  was  omitted.  Here  the  time  resolved  emissions  were  recorded 
and  used  as  the  background  for  the  experiments  of  interest.  It  was  observed  that  no 
background  corrections  were  necessary  for  Mixture  C (containing  0^)  whereas  Mixtures 
D and  E (containing  ^0)  needed  small  corrections.  The  calibration  factors  reported 
earlier  [14]  were  used  here.  Checks  made  during  the  course  of  these  experiments 
verified  those  values.  In  this  way  it  was  possible  in  effect  to  monitor  absolute 
oxygen  atom  concentrations.  All  data  were  collected  at  0.2-0. 5 ys  intervals,  and 

times  were  measured  relative  to  t . 

o 

The  results  of  a typical  experiment  with  Mixture  C is  shown  in  Fig.  2(a).  Four 
parameters  were  chosen  to  characterize  the  signal.  The  first  is  the  induction  time 
t.j  which  was  defined  as  the  time  at  which  the  signal  reached  0.02  V,  the  lowest  sig- 
nal which  could  be  reproducibly  extracted  from  the  data.  For  the  CO  concentration 

1 4 3 

used,  this  corresponds  to  an  oxygen  atom  concentration  equal  to  4.4x10  /cm  at 
14  3 

1705K  and  2.5x10  /cm  at  2455K.  After  the  induction  time  the  signal  rose  rapidly, 
achieved  a maximum  value,  and  then  decayed  at  a slower  rate.  The  initial  rise  may 
well  be  exponential,  but  it  does  not  continue  sufficiently  long  to  measure  an  un- 
equivocal growth  constant;  instead  the  time  t1  at  which  the  signal  reached  0.2  V 
Is  recorded.  The  maximum  voltage  reached  is  recorded  as  Vm,  and  the  time  at  which 

the  maximum  is  reached  is  listed  as  t . 

m 

Fig.  2(b)  Illustrates  the  results  observed  with  Mixtures  D and  E.  As  expected 
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from  the  earlier  work  [14],  no  induction  time  is  evident  here.  These  signals  are 
characterized  in  terms  of  two  parameters:  here  t1  designates  the  time  at  which  the 
signal  reached  one-half  its  maximum  value  and  Vm  is  the  maximum  signal.  (These 
parameters  are  also  listed  in  Table  2.)  The  flame-band  signal  here  does  not  exhibit 
the  obvious  decline  seen  in  Mixture  C,  and  itwasmore  difficult  to  assign  a t value 
here.  For  this  reason  no  t values  are  reported.  The  corresponding  plots  of  all  of 
the  parameters  in  Table  2 versus  temperature  exhibit  scatter  consistent  with  the  er- 
ror analysis.  (This  analysis  included  a 10%  uncertainty  in  the  calibration  factor.) 

Note  in  Table  2 that  the  temperature  dependence  of  Vm  is  much  greater  for  the 
N2O  systems.  This  same  phenomenon  was  observed  earlier  in  the  study  of  oxidation 
of  hydrogen  by  O2  and  N2O  [14]. 

CO2  Production 

The  Infrared  emissions  at  4.27  ym  were  collected  at  10.0  us  intervals.  Here 
all  time  data  were  measured  relative  to  the  window  midpoint  passage  time  (t  ). 
Notice  that  this  procedure  is  different  from  that  employed  for  the  flame-band  and 
CH2O  data.  Here  it  was  more  convenient  to  use  tmp  so  that  calculated  profiles  could 
be  directly  compared  to  those  observed.  Both  the  nature  of  the  CO2  signal  and  the 
slower  sampling  rate  here  meant  that  the  finite  width  of  the  observation  window  had 
little  effect;  no  explicit  corrections  need  be  applied. 

The  emission  observed  from  Mixture  C was  corrected  for  CO  emission  and  the 
resulting  signal  converted  to  an  absolute  CO2  concentration  as  in  the  earlier  work 
[14].  (More  recent  calibration  experiments  indicated  the  calibration  factor  used 
earlier  was  In  error  by  ^ 5%.)  A typical  C02  profile  is  shown  in  Fig.  3(a).  Here 
three  parameters  were  used  to  describe  the  signal:  (1)  The  induction  time  ti  was 
defined  as  the  Intercept  of  the  tangent  to  the  region  of  maximum  growth  r-?te  with 
the  time  axis.  (2)  The  maximum  growth  rate  was  labelled  S.  (3)  The  time  at  which 
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the  CO2  concentration  reached  4.5x10  /cm  was  listed  as  t".  The  parameters  are 
listed  in  Table  3. 

A typical  CO2  profile  for  an  ^0  experiment  is  shown  in  Fig.  3(b).  Here  an 

additional  correction  for^O  emission  must  be  made.  As  a result  it  was  felt  that 

15  3 

points  with  a C02  concentration  less  than  % 5x10  /cm  were  potentially  suspect. 

Thus  the  early  time  behavior  here  was  not  considered.  Instead  two  later  points 

15  3 

on  the  profile,  t’,  corresponding  to  the  time  when  [C02]  = 5x10  /cm  , and  t",  when 
[CO2]  = 1.5x10  /cm  were  used.  The  minimum  exhibited  for  t'1  in  both  D and  E at  in- 
termediate temperatures  is  similar  to  the  behavior  seen  earlier  in  the  H2  studies  [14]. 
The  maximum  CO2  production  was  seen  at  the  lower  temperatures  in  Mixtures  D and  E, 
whereas  in  C the  largest  CO2  production  came  at  the  highest  temperatures.  The  C02 
production  was  generally  lower  in  the  N2O  system;  this  is  reflected  in  the  different 
defining  relations  for  t".  As  with  the  other  reaction  parameters  discussed  above, 
an  error  analysis  of  the  CO2  data  indicates  the  observed  scatter  is  consistent  with 
the  assigned  errors.  In  some  cases,  the  error  bars  become  very  large;  after  the 
"knee"  in  the  CO2  profiles  small  calibration  errors  translate  into  large  time  errors. 


DISCUSSION 


As  outlined  In  the  Introduction,  the  purpose  of  this  study  was  to  attempt  to 
more  accurately  assign  high  temperature  rate  constants  to  the  reactions  of  formalde- 
hyde. The  approach  used  was  to  collect  a variety  of  data  on  CH20  systems  and  couple 
this  with  information  obtained  from  earlier  studies  on  simpler  systems  to  reduce  the 
number  of  variables  to  the  point  where  meaningful  assignments  could  be  made.  A 
plausible  mechanism  for  formaldehyde  oxidation  where  O2  is  the  oxidant  is  listed 
as  Reactions  (1)  through  (15)  in  Table  4.  Fortunately  quite  a bit  is  known  about 
many  of  these  reactions.  For  example,  Reactions  (9)  through  (15),  with  the  rate  con- 
stants listed,  have  been  shown  to  be  quantitatively  consistent  with  recent  experiments 
on  the  H2/02/C0  system  [14].  Likewise  CH20  pyrolysis  studies  [6]  have  yielded  appre- 
ciable information  about  Reactions  (1),  (2),  (5),  and  (6).  The  calculations  done 
during  this  work  have  shown  that  the  results  are  insensitive  to  the  values  used  for 
ky  and  kg.  Thus  the  goal  of  this  study  could  be  interpreted  as  trying  to  find  values 
of  kg  and  k4  that  would  satisfactorily  explain  the  reported  data  in  CH20/02  systems 
with  the  added  constraint  that  the  choices  for  k^,  k2,  kg,  and  kg  be  consistent  with 
the  earlier  pyrolysis  work. 

In  the  process  of  comparing  calculated  and  observed  profiles,  adequate  pro- 
vision must  be  made  for  both  the  finite  observation  window  width  and  the  finite  in- 
strument response  time.  The  results  described  below  were  obtained  by  the  following 
procedure:  (1)  The  system  of  rate  equations  was  first  numerically  integrated  in 
the  usual  way  [14]  to  obtain  concentration-time  profiles.  (2)  For  both  CH20  and 
flame-band  comparisons,  these  profiles  ([CH20]  and  [C0][0]  versus  time)  were  then 
integrated  over  a trapezoidal  slit  function  whose  dimensions  were  fixed  by  applica- 
tion of  simple  geometrical  optics  to  the  colllmatlon  system  used.  These  results,  in 
turn,  were  modified  to  account  for  the  finite  detector  response  time.  Here  the  detector 
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was  treated  as  the  capacitor  in  a simple  RC  circuit.  The  equation  describing  the 
capacitor  voltage  was  then  integrated  numerically  to  yield  a result  which  could  be 
compared  to  the  observations.  In  practice  it  was  observed  that  the  detector  response 
> could  be  adequately  treated  simply  by  adding  the  approximate  time  constant  (1  us  for 

IR,  2 ys  for  flame-band)  to  the  slit  routine  output.  The  modification  to  account 
for  the  finite  window  width  was  substantially  greater  than  the  subsequent  response 
time  correction.  In  the  flame-band  case,  the  calibration  factor  was  included  in  the 
calculations  so  the  final  output  was  simply  a voltage-time  profile  to  be  directly 
compared  to  that  observed;  there  was  no  need  for  any  normalization  here.  Given  the 
fact  that  the  CO  concentration  is  reasonably  constant,  this  comparison  in  effect  re- 
lates calculated  and  observed  oxygen  atom  concentrations  on  an  absolute  basis.  For 
CH2O,  there  is  no  such  calibration.  Here  the  observed  and  calculated  results  are 
normalized  at  the  plateau  (or  peak)  value.  This  method  of  comparison  of  slit-cor- 
rected output  to  the  data  dictated  the  choice  of  t as  the  reference  time  in  Tables 
1 and  2. 

For  CO2  comparisons,  slit  effects  are  of  no  consequence  since  the  first  meaning- 
ful data  comes  at  sufficiently  long  times  that  use  of  t is  adequate.  Here  the  cal- 
culated and  observed  profiles  could  be  directly  compared;  again  there  is  no  normali- 
zation required  since  absolute  calibration  factors  have  been  measured.  The  compari- 
sons were  normally  made  at  three  temperatures.  The  experimental  parameters  of  these 
temperatures,  with  associated  error  bars,  are  listed  in  Table  5.  In  most  cases  the 
values  of  the  various  parameters  have  been  taken  directly  from  the  appropriate  Table, 
but  in  a few  instances  the  values  listed  were  adjusted  slightly  to  better  conform  to 
the  other  data  in  that  temperature  range.  (In  cases  where  adjustments  were  made,  the 
amount  of  the  change  was  always  less  than  the  error  bar.) 

The  observations  were  first  compared  to  calculated  values  resulting  from  a choice 
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of  CH20  rate  constants  typical  of  those  presently  used  in  modeling  studies.  Set  I 
has  k]  = 8.3xl0'8  exp(-301  kJ/RT),  k2  = 3.3xl0-11  exp(-13.8  kJ/RT),  k3  = 8.3xl0'9 
exp(-54.3  kJ/RT),  and  k^  * 8.3x10"^  exp(-19.2  kJ/RT).  These  values  were  used  in 
the  recent  study  of  Olson  and  Gardiner  [4].  Rate  constants  for  Reactions  (5)  - (15) 
were  those  listed  in  Table  4.  Note  this  Set  (Table  5)  yields  calculated  values  in 
marked  disagreement  with  those  observed.  Both  CH20  decay  in  A and  oxygen  atom  pro- 
duction in  C is  much  too  fast.  Vm  at  2455K  is  also  much  too  high.  Likewise  the  C02 
profiles  In  C don't  fit  well.  Since  the  recent  pyrolysis  study  [6]  indicated  that 
k1  was  much  lower  than  the  value  used  in  Set  I,  a series  of  calculations  were  done 
using  the  k^  and  k2  values  from  Ref.  [6]  in  conjunction  with  the  k3  and  k4  values 
of  Set  I.  Although  some  improvement  was  seen,  these  results  suggested  the  rates  were 
still  larger  than  observed.  These  results  suggested  that  the  k3  and/or  k4  values 
were  too  high  in  Set  I. 

It  is  likely  that  the  lower  limits  of  k3  and  k^  at  high  temperatures  are  those 
represented  by  a simple  Arrhenius  extrapolation  of  the  low  temperature  data.  These 
values  (listed  in  Table  4)  have  both  been  reported  very  recently,  and  another  calcula- 
tion (Set  II)  was  done  with  these  values  and  the  k^  and  k2  values  of  Ref.  6 [k-j  = 

6x10  ^ exp(-364  kJ/RT)  and  k2  = 4.2x10”^  exp(-16.0  kJ/RT)]  which  fit  the  pyrolysis 
data.  Here  there  was  further  improvement,  but  there  were  still  some  problems  at  the 
high  temperatures  in  Mixture  C where  this  combination  still  appeared  to  be  too  fast 
(note  In  particular  that  t'-t.  for  the  flame-band  was  only  4 ps  vs.  7 observed  and 
that  the  maximum  rate  of  C02  production  was  higher  than  observed).  This  conclusion 
was  verified  by  explicit  comparison  of  complete  concentration-time  profiles. 

The  earlier  pyrolysis  study  Indicated  that  one  could  achieve  a satisfactory 
fit  there  with  a variety  of  k1 , k2  choices  as  long  as  the  product  of  the  values  re- 
mained approximately  constant.  Additional  calculations  suggested  that  the  rate  constant 
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combination  designated  in  Table  5 as  Set  III  (k1  = 5.5xl0-8  exp(-339  kJ/RT),  k2=5.5xl0‘10 
exp(-43.9  kJ/RT),  k^  and  k^  same  as  Set  II)  gave  a reasonable  fit  to  both  the  CH^O 
oxidation  and  pyrolysis  studies.  Note  in  Table  5 that  most  of  the  calculated  para- 
meters now  lie  within  the  experimental  error  bars.  Again,  explicit  profile  compari- 
sons confirm  the  adequacy  of  the  fit.  Fig.  1 illustrates  this  fit  for  Mixture  A and 
Figs.  2(a)  and  3(a)  demonstrate  the  fit  for  Mixture  C.  Of  the  three  parameter  values 
that  don't  fall  within  the  error  bars,  there  appear  to  be  reasonable  explanations  for 
the  deviations.  The  most  serious  deviation  is  seen  for  Mixture  A at  1670  where  the 
calculated  t"  is  too  early.  A possible  explanation  for  this  discrepancy  was  dis- 
cussed earlier;  t"  measurements  may  be  unduly  influenced  by  the  background  emissions 
here.  Note  the  low  temperature  data  on  Mixture  C are  in  reasonable  agreement.  It 
would  appear  that  the  only  way  one  could  slow  down  the  calculation  to  fit  the  reported 
t"  value  at  1670  would  be  either  to  use  k^  and/or  k^  values  lower  than  used  here  or 
use  a k^,  k2  combination  which  would  not  be  consistent  with  the  pyrolysis  results. 

Given  the  definitive  nature  of  the  low  temperature  studies,  it  is  unlikely  that  such 
an  adjustment  in  k^  or  k^  is  justified.  It  would  appear  then  that  the  fault  here  lies 
with  the  data.  In  a similar  light,  the  slightly  higher  maximum  rates  calculated  at 
1705  and  1940K  appear  to  present  no  real  problem.  Here  the  calculated  parameters  are 
just  slightly  out  of  the  error  bars  and  the  measurements  are  difficult  since  data  were 
only  collected  at  10  ys  intervals. 

The  rate  constant  combinations  used  in  Set  III  would  thus  appear  to  be  remarkably 
successful  at  explaining  a large  body  of  disparate  data  on  CH20  systems.  (The  k-j , 
k2  values  of  Set  III  also  quantitatively  describe  the  pyrolysis  results  of  Ref.  6j 
A particularly  satisfying  aspect  of  the  assignments  is  that  very  few  rate  constants 
needed  adjustment.  In  essence  the  only  variables  were  k^  and  k2-  Note  k^  and  k^ 
were  taken  directly  from  low  temperature  data.  There  is  still  some  uncertainty  with 
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respect  to  kg  and  kg.  Recently  Reilly,  et  al . [24]  suggested  a room  temperature 
value  of  kg  a factor  of  two  larger  than  that  used  here.  The  pyrolysis  studies  [6] 

Indicated  simply  that  the  ratio  used  in  that  work  was  needed  to  achieve  a reasonable 
fit.  If  the  kg  value  used  here  is  indeed  somewhat  low,  this  would  simply  indicate 
kg  is  also  somewhat  low;  other  conclusions  would  remain  unaffected. 

A sequence  of  calculations  using  various  k^ , k2  values  on  Mixture  C suggests 
that  one  can  swing  from  one  end  of  the  error  bars  to  the  other  with  less  than  a fac- 
tor of  two  changes  in  these  rate  constants.  It  would  appear  then  that  the  values 
specified  in  Set  III  are  probably  accurate  to  within  a factor  of  two  within  the  1700- 
2500K  temperature  range.  Similarly,  studies  on  Mixture  A suggest  that  the  k^,  k^ 
values  used  here  are  probably  accurate  to  within  a factor  of  two  over  this  interval. 

An  additional  check  upon  these  rate  constant  assignments  could  be  obtained  from  \ 

analysis  of  the  oxidation  of  C^O  by  N20.  Here  the  concentration-time  behavior  of  1 

the  various  reactive  species  is  radically  different.  Another  advantage  of  using  this 
system  for  analysis  is  that  no  additional  unknowns  are  generated.  The  ^0  dissocia- 
tion kinetics  are  reasonably  well  characterized  [20,25]  and  its  use  as  an  oxidant 
for  Hg  has  also  been  studied  [14,26]. 

Earlier  it  was  pointed  out  that  a small  change  was  made  in  the  C02  calibration 
factor  for  this  work.  As  a result,  the  earlier  hydrogen  work  [14]  was  reexamined 
to  see  if  any  rate  constant  assignments  made  then  needed  to  be  changed.  An  examina- 
tion of  the  H2/O2/CO  system  Indicated  no  changes  were  mandated.  A better  description 
of  the  Hg/NgO/CO  system  could  now  be  obtained  with  a slightly  larger  value  of  the 
rate  constant  for  H + N20  = N2  + OH.  This  new  value  is  klg  = 1 . 5x1 0”9  exp(-92  kJ/RT); 
this  value  compares  favorably  with  that  reported  at  lower  temperatures  by  Albers, 
et  al.  [27].  There  is  still  evidence  of  non-Arrhenius  behavior,  but  not  as  much  as 
previously  reported.  This  value  of  k^g  was  Incorporated  into  the  present  analysis. 
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The  other  NgO  rate  constants  (Table  4)  were  taken  from  the  earlier  work  [14], 

The  N2O  studies  were  conducted  analogously  to  the  work.  In  Mixture  B, 

C^O  decay  was  monitored,  and  CO2  and  oxygen  atom  production  were  followed  in  Mix- 
tures D and  E.  Selected  experimental  parameters  are  listed  in  Table  6 and  compared 
to  various  calculated  results.  Here  Set  I of  the  calculations  used  the  same  values 
used  in  Set  I of  the  02  system.  Again  some  mismatches  are  evident,  but  the  differences 
tend  not  to  be  as  dramatic  as  those  seen  with  the  02  system.  Use  of  Set  III  (the 
set  which  worked  well  for  the  O2  system)  again  yields  encouraging  results.  Good 
agreement  is  seen  with  Mixture  B.  In  Mixtures  D and  E,  two  problems  appear:  (1)  The 
calculated  Vm  is  somewhat  high  at  higher  temperatures.  (2)  CO2  production  seems  a 
little  slow  (particularly  at  the  lower  temperatures).  It  is  interesting  to  note  that 
in  each  case  Mixture  E presents  more  of  a problem.  The  causes  of  these  discrepancies 
are  not  understood.  It  is  possible  that  they  are  experimental  in  origin.  On  the  one 
hand,  the  C^O  concentrations  were  very  low  here  and  there  may  well  be  an  adsorption 
and/or  desorption  problem  --  changes  induced  by  these  events  would  have  a much  greater 
effect  in  these  low  concentration  mixtures.  An  alternative  explanation  is  that  there 
are  more  difficult  background  corrections  to  be  made  in  the  N20  system.  Of  course, 
there  is  always  the  possibility  of  some  unexpected  mechanistic  complexities.  In 
spite  of  some  problems  here,  the  overall  agreement  is  still  quite  reasonable.  The 
explicit  profile  comparisons  shown  in  Figs.  2(b)  and  3(b)  reinforce  this  point.  Par- 
ticularly noteworthy  is  that  the  fit  achieved  was  obtained  with  absolutely  no  varia- 
tion in  any  rate  constants.  This  fact  seems  to  reinforce  the  basic  accuracy  of  the 
assignments  made  in  the  oxygen  system. 

Reactions  (1)  - (15)  with  the  rate  constants  listed  in  Table  4 were  then  used 
In  a series  of  calculations  simulating  the  experiments  of  Gay,  et  al . [5],  As  expected 
from  the  results  of  Ref.  6,  the  calculated  results  for  the  pyrolysis  experiments  sug- 


gested  a stronger  temperature  dependence  than  reported;  reasonable  agreement  was 
seen  at  1950K,  but  the  calculated  results  were  much  slower  at  1550K.  Conversely, 
in  the  oxidation  experiments,  there  was  reasonable  agreement  at  1550K  with  the  cal- 
culated results  being  too  fast  at  1980K.  One  interesting  point  is  that  the  calcu- 
lated CH2O,  CO,  Hg,  and  ^0  profiles  at  1550K  are  in  surprisingly  good  agreement 
with  those  shown  in  Fig.  2 in  Ref.  5.  Even  at  1 71  OK,  the  reported  time  for  20% 
reaction  of  C^O  of  120  ys  compares  favorably  with  the  calculated  value  of  112  us. 

By  1980K,  however,  the  calculated  time  is  only  32  ys  as  opposed  to  the  reported  value 
of  70  ys.  In  spite  of  differences  in  temperature  dependencies  it  would  appear  that 
the  present  mechanism/rate  constant  combination  predicts  results  reasonably  consistent 
with  the  oxidation  studies  of  Gay,  et  al.  It  is  significant  to  note  that  the  experi- 
mental conditions  of  that  work  differed  significantly  from  that  reported  here  ( i . e . , 
lower  temperatures,  lower  total  concentrations,  higher  reactant  concentrations,  and 
different  diagnostics). 
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CONCLUSIONS 


The  formaldehyde  oxidation  results  presented  in  this  work  are  interpreted  in 
terms  of  formaldehyde  reaction  rate  constants  appreciably  lower  than  usually  seen 
in  the  high  temperature  literature.  For  example,  formaldehyde  dissociation  here 
Is  reported  to  be  an  order  of  magnitude  slower  than  reported  by  Schecker  and  Jost  [2] 
and  several  orders  of  magnitude  lower  than  Peeters  and  Mahnen  [3]  suggest.  Similarly, 
the  rate  of  hydroxyl  attack  on  formaldehyde  is  at  least  an  order  of  magnitude  lower 
than  values  currently  used  in  some  modeling  studies  [4,9].  The  values  for  hydrogen 
atom  and  oxygen  atom  attack  are  much  closer  to  those  now  used;  the  oxygen  rate  con- 
stant is  lower  by  a factor  of  two  while  the  hydrogen  atom  value  is  higher  by  about 
this  same  amount.  One  interesting  aspect  of  the  results  is  that  0 and  OH  reactions  1 

appear  to  be  best  interpreted  in  terms  of  a simple  Arrhenius  extrapolation  of  lower 
temperature  data  while  the  hydrogen  atom  reaction  appears  to  have  some  non-Arrhenius 
character  when  compared  to  the  recent  results  of  Klemm  [28].  | 

One  possible  reason  for  the  dramatic  differences  between  this  work  and  the  liter- 
ature could  well  be  due  to  a lack  of  understanding  of  the  methane  oxidation  mechanism. 

Until  the  recent  suggestions  [9,29]  that  C2  species  may  play  a significant  role  in 
methane  oxidation,  most  mechanistic  studies  suggested  the  sequence  CH4  -*■  CH3  -*•  CH20  -*• 

HCO  -►  CO  -*■  C02.  Given  the  lack  of  direct  high  temperature  information  about  formal- 
dehyde reactions,  explanations  based  upon  decay  of  CH^  or  formation  of  CO  or  C02  could 
easily  assign  too  large  a rate  to  CH20  reactions  rf  an  alternative  path  were  actually 
present  but  not  accounted  for  in  the  analysis.  The  recent  studies  suggesting  such 
analternate  path  may  well  be  consistent  with  the  present  study  which  could  simply 
Indicate  (In  this  context)  that  all  of  the  CH4  doesn't  follow  the  CH20  sequence. 

In  this  light,  the  present  work  highlights  the  dangers  inherent  in  mechanistic  analysis 


in  complex  systems  where  there  are  simply  too  many  variables  for  adequate  character- 
ization. 
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TABLE  1.  SUMMARY  OF  OBSERVATIONS  ON  CH20  DECAY 


Mixture3 

T/K 

M/1018cm-3 

iiMc 

t'/psd 

t"/y 

2455 

4.62 

• 

11 

17 

A 

2400 

4.63 

- 

10 

17 

2370 

4.55 

- 

13 

20 

0.50%CH,0 

2235 

4.49 

- 

13 

22 

C 

2185 

4.37 

- 

16 

26 

0.88%0, 

2120 

4.36 

- 

21 

33 

c 

2005 

4.31 

23 

30 

50 

1980 

4.14 

22 

31 

53 

1970 

4.24 

22 

32 

55 

1935 

4.09 

28 

36 

64 

1850 

4.14 

37 

45 

84 

1815 

4.06 

46 

50 

102 

1780 

4.09 

47 

50 

110 

1750 

4.06 

60 

72 

136 

1670 

4.09 

78 

- 

175 

B 

2275 

4.73 

11 

19 

2225 

4.70 

- 

10 

20 

1.01XCH,0 

2165 

4.74 

- 

12 

19 

2130 

4.65 

- 

13 

21 

1.07%N,0 

2120 

4.61 

- 

13 

25 

C 

2040 

4.48 

- 

17 

32 

1990 

4.40 

- 

23 

45 

1965 

4.36 

- 

20 

42 

1940 

4.42 

- 

18 

40 

1900 

4.39 

19 

25 

65 

1895 

4.38 

18 

24 

55 

1785 

4.26 

22 

32 

82 

1755 

4.27 

26 

36 

96 

1740 

4.22 

30 

34 

124 

1685 

4.09 

40 

40 

148 

aBalance  of  mixture  was  argon, 
kjotal  concentration 

cT1me  relative  to  t^,  (see  text)  for  signal  to  decay  5*  from  an  early  constant  signal. 

dFor  Mixture  A,  time  relative  to  t0  when  the  CH2O  emission  drops  to  3.0  volts.  For  B, 
time  when  the  emission  drops  to  1.0  volts. 

eFor  A,  time  when  emission  drops  to  0.3  V;  for  B,  time  when  emission  drops  to  0.1  V. 


TABLE  2.  SUMMARY  OF  FLAME-BAND  OBSERVATIONS 


Mixture3 

T/K 

M/1018cm'J 

iiZH£C 

t'/psd 

tmM6 

Vm/ volts 

C 

2455 

5.14 

6 

13 

63 

1 .63 

2380 

5.15 

5 

13 

70 

1.55 

0.50%CH,0 

2290 

5.08 

8 

19 

78 

1.33 

c 

2265 

5.07 

6 

17 

78 

1.40 

0. 92%0o 

2160 

4.93 

12 

23 

85 

1.15 

c 

2025 

4.93 

15 

37 

96 

1.08 

9.1 5%C0 

1950 

4.74 

27 

46 

124 

1.03 

1940 

4.64 

32 

54 

140 

0.90 

1940 

4.58 

33 

55 

135 

0.95 

1845 

4.55 

61 

87 

172 

0.78 

1810 

4.56 

74 

102 

200 

0.80 

1800 

4.51 

70 

103 

180 

0.77 

1735 

4.42 

71 

104 

190 

0.75 

1705 

4.29 

82 

126 

205 

0.67 

1660 

4.28 

148 

190 

285 

0.60 

1605 

4.32 

219 

275 

385 

0.60 

D 

3010 

5.59 

5 

2.51 

2890 

5.49 

7 

2.10 

0.049%CHo0 

2825 

5.42 

7 

2.15 

c 

2805 

5.45 

5 

2.08 

1.02%No0 

2680 

5.33 

7 

1.64 

c 

2540 

5.25 

10 

1.33 

8.21 %C0 

2440 

5.15 

12 

1.17 

2315 

5.08 

17 

0.87 

2215 

4.99 

21 

0.63 

2180 

4.90 

. 25 

0.65 

2125 

4.87 

24 

0.49 

2060 

4.77 

37 

0.37 

2010 

4.74 

52 

0.30 

1950 

4.66 

52 

0.22 

1905 

4.78 

63 

0.17 

(continued) 
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TABLE  2.  SUMMARY  OF  FLAME-BAND  OBSERVATIONS  (continued) 


Mixture* 

I/k 

M/1018cnf3 

iiMc 

f/psd 

.Wh?.6 

V^volts1 

E 

3045 

5.54 

5 

2.64 

2955 

5.49 

6 

2.55 

0.026%CH90 

2895 

5.47 

7 

2.43 

L 

2850 

5.41 

6 

2.26 

1 .02%N90 

2705 

5.36 

8 

1.92 

c 

2650 

5.28 

9 

1.77 

8.10%C0 

2575 

5.25 

11 

1.68 

2470 

5.15 

15 

1.33 

2385 

5.06 

15 

1.17 

2145 

4.88 

30 

0.72 

2085 

4.82 

39 

0.48 

1940 

4.66 

74 

0.30 

flBalance  was  argon. 
bTotal  concentration 

cTime  relative  to  tn  (see  text)  when  signal  * 0.02  A/. 

j o 

For  C,  time  when  signal  = 0.20  volts.  For  D and  E,  time  when  signal  is  one-half 
maximum  value. 

eT1me  relative  to  tQ  when  maximum  signal  achieved. 

^Maximum  signal. 
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TABLE  3.  SUMMARY  OF  C02  OBSERVATIONS 


Mixture® 

T/K 

M/1018cm“3 

t c 
— i— 

S/1020cm"3s_1 

t'/use 

t'Vus 

C 

2455 

5.14 

6 

19.5 

32 

2380 

5.15 

8 

19.8 

35 

0. 50%CH,0 

2290 

5.08 

11 

16.0 

44 

2265 

5.07 

8 

15.7 

41 

0.92%0, 

2160 

4.93 

15 

12.6 

58 

2140 

5.01 

14 

12.7 

56 

9.1 5%C0 

2025 

4.93 

21 

9.26 

79 

1950 

4.74 

32 

8.68 

97 

1940 

4.64 

34 

6.79 

120 

1940 

4.58 

37 

6.66 

120 

1845 

4.55 

60 

4.61 

200 

1810 

4.56 

75 

4.67 

202 

1800 

4.51 

76 

4.49 

208 

1735 

4.42 

74 

4.32 

209 

1705 

4.29 

92 

4.19 

238 

1660 

4.28 

146 

2.83 

370 

1605 

4.32 

218 

2.32 

560 

D 

3010 

5.59 

5 

325 

2890 

5.49 

8 

235 

0.049%CHo0 

2825 

5.42 

5 

160 

c 

2805 

5.45 

5 

112 

1 .022N90 

2680 

5.33 

5 

65 

c 

2540 

5.25 

7 

25 

8.2UC0 

2440 

5.15 

8 

23 

2315 

5.08 

9 

25 

2215 

4.99 

13 

33 

2180 

4.90 

13 

32 

2125 

4.87 

14 

37 

2060 

4.77 

22 

54 

2010 

4.74 

25 

63 

1.950 

4.66 

36 

82 

1905 

4.78 

45 

106 

(continued) 
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TABLE  3.  SUMMARY  OF  C02  OBSERVATIONS  (continued) 


Mixture* 

m 

M/1018cm'3b 

t^  S/1020cnfVld 

t !Me 

t'Vu 

E 

3045 

5.54 

13 

557 

2955 

5.49 

10 

475 

0.026%CH,0 

2895 

5.47 

10 

425 

c 

2850 

5.41 

9 

380 

1.02%N,0 

2705 

5.36 

8 

278 

2650 

5.28 

. 10 

278 

8.10%C0 

2575 

5.25 

9 

200 

2470 

5.15 

10 

60 

2385 

5.06 

11 

47 

2260 

4.99 

15 

46 

2145 

4.88 

20 

52 

2085 

4.82 

27 

72 

1940 

4.66 

42 

105 

aBalance  was  argon. 

^Total  concentration. 

cTime  relative  to  t (see  text)  when  tangent  to  [C0,]-time  plot  intercepts  t-axis. 
Slope  of  initial  region  (after  induction  period)  of  plot  of  [C02]  vs.  time. 
eTime  relative  to  tm_  when  [C0o]  = 5xl0^/cm3. 

f mP  t 1 c o 

TFor  Mixture  C,  time  relative  to  t_  when  [CO,]  s 4.5x10  /cm  ; for  D and  E,  time 
when  [C02]  » 1.5xl016/cm3.  mp  L 


TABLE  4.  MECHANISM  AND  RATE  CONSTANTS 


Rate  Constant3 


Reaction 

1^3.10-^ 

ii 

h 

Reference 

1. 

CH20  + M 

= 

HCO  + H + M 

-7.26 

0 

339 

this  work 

2. 

CH20  + H 

= 

HCO  + H2 

-9.26 

0 

44 

this  work 

3. 

CH20  + OH 

= 

HCO  + H20 

-10.90 

0 

OJ 

15 

4. 

ch2o  + 0 

= 

HCO  + OH 

-10.52 

0 

12.9 

16 

5. 

HC0  + M 

= 

H + CO  + M 

-9.59 

0 

61.4 

2 

6. 

HCO  + H 

= 

h2  + CO 

-9.48 

0 

0 

7 

7. 

HCO  + OH 

= 

h2o  + CO 

-9.77 

0 

0 

7 

8. 

HCO  + 0 

= 

OH  + CO 

-9.77 

0 

0 

7 

9. 

H + 02 

= 

OH  + 0 

-6.69 

-.91 

69.5 

17 

10. 

0 + h2 

= 

OH  + H 

-9.44 

0 

57.5 

18 

11. 

o + h2o 

= 

OH  + OH 

-9.96 

0 

76.8 

19 

12. 

0 + CO  + M 

= 

C02  + M 

-33.80 

0 

0 

20 

13. 

OH  + H2 

= 

H20  + H 

-10.44 

0 

21.5 

21 

14. 

OH  + CO 

= 

C02  + H 

-11.18 

0 

33.4 

22 

15. 

CO  + o2 

= 

co2  + 0 

-11.24 

0 

209 

23 

16. 

N20  + M 

= 

N2  + 0 + M 

-9.57 

0 

216 

14 

17. 

n2o  + 0 

s 

NO  + NO 

-10.11 

0 

117 

20 

18. 

n2o  + 0 

— 

N?  + 07 

-10.11 

0 

117 

20 

19. 

A 

N20  + H 

OH  + N2 

-8.82 

0 

1 

92 

this  work 

Expressed  in  the  form  k * ATnexp(-E^/RT)  in  cm3,  molecule,  sec,  and  kJ  units. 
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TABLE  5.  COMPARISON  OF  CALCULATIONS  AND  OBSERVATIONS  (02  SYSTEMS) 


Mixture3  Parameter1* 

T/K 

Observed 

Set  I c 

Set  11^ 

Set  III 

A ti 

1670 

78±11 

32 

66 

68 

1935 

28±  6 

- 

23 

24 

t ' 

1935 

36±  5 

a. 

34 

36 

2455 

11±  2 

10 

11 

11 

t" 

1670 

175±21 

74 

137 

140 

1935 

64±12 

- 

65 

66 

2455 

17±  2 

13 

18 

19 

C ti(F.B.) 

1705 

97±24 

46 

78 

86 

1940 

32±  5 

- 

28 

34 

2455 

6±  2 

5 

6 

8 

t'(F.B.) 

1705 

135±20 

79 

117 

125 

1940 

54±  6 

- 

49 

56 

2455 

13±  2 

8 

10 

14 

tm(F.B.) 

1705 

228±28 

180 

212 

222 

1940 

135±21 

- 

112 

122 

2455 

63±20 

48 

48 

51 

Vm(F.B.) 

1705 

0.67±.09 

0.61 

0.59 

0.61 

1940 

0.94±.10 

- 

0.89 

0.90 

2455 

1 . 63± .17 

2.25 

1.85 

1.65 

t^COg) 

1705 

100±13 

53 

83 

94 

1940 

35±  5 

- 

33 

40 

• 

2455 

6±  3 

4 

5 

8 

S(C02) 

1705 

3.7±0.6 

3.8 

4.2 

4.4 

c 

1940 

7.4+0. 9 

- 

8.5 

8.8 

2455 

20 . 6±3 . 5 

22.6 

25.8 

22.7 

t"(co2) 

1705 

270±32 

209 

229 

239 

c 

1940 

1 1 5±1 5 

• 

102 

107 

2455 

32±  4 

47 

33 

35 

aSee  Tables  1 and  2 for  composition. 

bSee  Tables  1,  2 and  3 for  definitions. 


c" Literature  values";  see  text  for  details. 
d"Or1g1nal  Pyrolysis  Fit";  see  text  for  details. 
e"Mod1f1ed  Pyrolysis  Fit";  see  text  for  details. 


TABLE  6.  COMPARISON  OF  CALCULATIONS  AND  OBSERVATIONS  (N20  SYSTEMS) 


Mixture3 

Parameter15 

I/K 

B 

*1 

1740 

1900 

f 

1740 

1900 

2275 

t" 

1740 

1900 

2275 

D 

t' (F.B. ) 

1950 

2540 

2890 

VJF.B.) 

1950 

m 

2540 

2890 

t'(co2) 

1950 

2540 

2890 

t"(C02) 

1950 

2540 

2890 

E 

t1 (F.B. ) 

1940 

2470 

2955 

V (F.B.) 
m 

1940 

2470 

2955 

t'(co2) 

1940 

2470 

2955 

t"(C02) 

1940 

2450 

2955 

a$ee  Tables  1 

and  2 for  composition. 

“Defined  In  Tables  1,  2,  and  3. 

c"L1terature  values";  see  text. 

"Modified  Pyrolysis  Pit";  see  text. 

Observed 

Set  I c 

Set  md 

30±  7 

25 

25 

19±  3 

- 

17 

36±  5 

29 

31 

24+  3 

- 

23 

11±  2 

10 

11 

115+21 

93 

109 

60±11 

- 

54 

19±  3 

13 

19 

56±  8 

47 

47 

10±  2 

- 

12 

7+  2 

9 

9 

0.22±.04 

0.26 

0.23 

1.37±.17 

- 

1.47 

2.24±.25 

2.53 

2.57 

36±  7 

45 

41 

7±  3 

- 

10 

6+  3 

5 

6 

82±13 

109 

84 

25±12 

- 

29 

21 5±55 

197 

221 

74±13 

60 

61 

14±  3 

- 

15 

6+  3 

9 

9 

0.30+.05 

0.33 

0.30 

1.33+.14 

- 

1.51 

2.55±.29 

3.03 

3.06 

42±  4 

60 

54 

10+  3 

- 

16 

10+  5 

9 

11 

1 05±1 0 

159 

125 

80±29 

- 

124 

475±63 

625 

641 

32 


LEGENDS 

Fig.  1.  A typical  CHgO  decay  plot  of  the  log  (corrected  I.R.  signal)  vs.  time. 

Here  tQ  is  the  time  at  which  the  reflected  shock  front  first  enters  the 
detector's  field  of  view.  The  open  circles  designate  every  fourth  data 
point  obtained  with  Mixture  A at  1935K;  the  closed  circles  are  the  cal- 
culated results  using  Set  III.  (See  text.) 

Fig.  2.  Typical  flame-band  emission  profiles.  Here  tQ  is  the  time  at  which  the 
reflected  shock  front  first  enters  the  detector's  field  of  view.  The 
open  circles  designate  every  twentieth  data  point  observed;  the  closed 
circles  represent  calculated  profiles  using  Set  III.  (See  text.) 

(a)  Mixture  C,  T 3 1940K,  M 3 4.64xl018/cm3. 

(b)  Mixture  D,  T 3 1950K,  M 3 4.66xl018/cm3. 

Fig.  3.  Typical  C02  production  profiles.  Here  t is  the  time  the  reflected 

shock  front  passes  the  center  of  the  observation  window.  The  open  cir- 
cles designate  experimental  data;  closed  circles  represent  calculated 
profiles  using  Set  III.  (See  text.) 

(a)  Mixture  C,  T 3 1940K,  M 3 4.64xl0^8/cm3. 

(b)  Mixture  D,  T 3 1950K,  M 3 4.66xl018/cm3. 
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and  for  CO?  production  via  IR  emission  at  4.27  These  data  were  then  com- 
pared to  the  results  obtained  from  numerical  integration  of  the  rate  equations 
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the  formaldehyde  reaction  rate  constants  listed  below: 

(1)  CH?0  + M = HCO  + H + M,  kj  = 5.5x10""  exp(-339  kJ/RT)  cm3  s'1 

(2)  CH20  + H = HCO  + H2,  k2  = 5.5xlO'10  exp(-43.9  kJ/RT)  cm3  s"1 

(3)  CH20  +011  = HCO  + H?0,  k3  = 1 .25x10”  1 ' exp ( - 0 - 7 kJ/RT)  cm3  s"’ 

(4)  CH20  + 0 = HCO  + OH,  k4  = 3.0xl0'U  exp(-12.9  kJ/RT)  cm3  s"1 

The  values  listed  for  k3  and  ki,  are  taken  directly  from  the  low  temperature 
studies  of  these  reactions;  the  k\  and  k2  values  also  fit  recent  experiments  in 
this  laboratory  on  CH20  pyrolysis.  The  N20  experiments  also  served  as  a check 
on  the  rate  constant  assignment  for  the  reaction: 

(19)  H + N?0  = N2  + OH,  k]9  = 1.5x10"°  exp(-92  kJ/RT)  cm3  s"1. 

Reasonable  agreement  was  achieved  here.  Sensitivity  studies  would  suggest 
these  rate  constants  are  probably  accurate  to  within  a factor  of  two  between 
1 700-2500K. 

The  values  of  kj  and  k3  obtained  in  this  work  are  at  least  an  order  of 
magnitude  lower  than  the  values  frequently  used  in  current  modeling  studies  of 
methane  oxidation.  The  implications  of  these  differences  are  discussed  in 
light  of  recent  suggestions  that  C2  species  participate  in  methane  oxidation. 
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